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Abstract 

A scale invariant spectrum of isocurvature perturbations is generated during col- 
lapse in the scaling solution in models where two or more fields have steep negative 
exponential potentials. The scale invariance of the spectrum is realised by a tachyonic 
instability in the isocurvature field. We show that this instability is due to the fact 
that the scaling solution is a saddle point in the phase space. The late time attractor 
is identified with a single field dominated ekpyrotic collapse in which a steep blue spec- 
trum for isocurvature perturbations is found. Although quantum fluctuations do not 
necessarily to disrupt the classical solution, an additional preceding stage is required 
to establish classical homogeneity. 

1 Introduction 

Given the success of simple inflation models in explaining the primordial homogeneity and 
almost scale-invariant spectrum of curvature perturbations on large scales, it is important 
to consider whether there is any other model for the early universe that might provide an 
alternative explanation. One contender is the idea of an early collapse phase rather than 
an accelerated expansion. This invokes the same basic mechanism as inflation to generate 
large-scale perturbations. A shrinking comoving Hubble scale can take quantum vacuum 
fluctuations of scalar fields on small scales and evolve them into overdamped perturbations 
on super-Hubble scales. But the approach to solving the homogeneity problem is quite 
different. 

In pre big bang model of Gasperini and Veneziano pQ the collapsing universe (in the 
Einstein frame [2]) is dominated by the kinetic energy of massless fields in the low energy 
string effective action. Such a solution is an attractor in a collapsing universe for a wide 
range of interaction potentials so offers a solution to the homogeneity problem, though this 
kinetic-dominated collapse is only marginally stable with respect to anisotropic shear [3J. By 
contrast the ekpyrotic scenario proposed an ultra-stiff collapse phase dominated by a steep 
negative potential [I] (see also [5j[6]). This is a late time attractor and is stable with respect 
to anisotropy. 
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On the other hand both models produce a steep blue spectrum of perturbations (in the 
spatially flat gauge) in the fields driving the collapse and hence in the comoving curvature 
perturbation [7J [8j [DJ [10] . A collapse phase dominated by a field with a pressureless equation 
of state does yield a scale-invariant spectrum of comoving curvature perturbations [TTl [T21 [T3] 
but it is not a late time attractor [TH [T5] . 

It was shown in [T6j that axion fields in the pre big bang can acquire a scale-invariant 
spectrum of isocurvature field perturbations in a pre big bang phase and that this could 
in principle be converted to curvature perturbations at a bounce [T7] or some time after 
the bounce via the curvaton mechanism [18[ HI]- The idea of producing a scale-invariant 
spectrum of isocurvature perturbations in axion fields was also investigated in the ekpyrotic 
scenario [20| [21|. 

In the ekpyrotic scenario it is also possible to generate a scale-invariant spectrum of 
isocurvature perturbations in models where two or more fields have steep exponential po- 
tentials. In this case there exists a scaling solution [22j [23] where the energy densities of the 
fields grow at the same rate during collapse, analogous to assisted inflation solutions in an 
expanding universe [2l] . Several authors [251 (2S1 [27] have recently proposed new ekpyrotic 
models based on this scaling solution. 

If a scale- invariant spectrum of comoving curvature perturbations can be generated during 
the collapse phase then there is still the "graceful exit problem" of turning collapse to 
expansion. In the pre big bang scenario it was envisaged that such a bounce could be 
achieved by higher-order corrections to string effective action which violate the null energy 
condition [2SJ [2H [3D]. A bounce due to higher-order kinetic terms in a ghost condensate 
model has recently been considered in the new ekpyrotic scenarios [2S1E7J (see also [3Tj 132"]). 

In this paper we discuss the tachyonic instability of long wavelength perturbations about 
the new ekpyrotic scaling solution which leads to the scale- invariant spectrum [22 | [25l |26[ [27] . 
We show that the solution is not an attractor at early or late times, but rather a saddle point 
in the phase-space. The late time attractor is the old ekpyrotic collapse dominated by a single 
field and the instability drives the scaling solution to this late time attractor. By contrast 
we find steep blue spectra for isocurvature field perturbations about the kinetic dominated 
early time attractors, or single-field dominated ekpyrotic late time attractors. We present 
the phase space analysis in an Appendix. 



2 Homogeneous field dynamics 

During the ekpyrotic collapse the expansion of the universe is assumed to be described by 
a 4D Friedmann equation in the Einstein frame with scalar fields with negative exponential 
potentials 

3H 2 = V + ^l (1) 

where 

V = - £ Vie-** , (2) 

i 

and we take Vi > and set 87rG equal to unity. 
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The authors of [25] found a scaling solution (previously studied in (22j [23]) in which both 
fields roll down their potential as the universe approaches a big crunch singularity, analogous 
to the assisted inflation dynamics found by [23] for positive potentials in an expanding 
universe. In this assisted ekpyrotic collapse we find power-law solution for the scale factor 

2 

a oc (— t) p , where p = — , (3) 

t C i 

where 

rh 2 —V-p-^fo C 2 

<P 2 -Vje-^i cf ' y ! 

As in the case of assisted inflation, this behaviour is easily understood after a rotation 
in field space [221 E3] into so-called adiabatic field, a, along the background trajectory (j4]), 
and the orthogonal, entropy field directions, Xii [33] . The potential ([2]) can then be simply 
re-written as 

V = -U( Xi ) e-~ , (5) 

where 

° i H 

To simplify the analysis we restrict ourselves to the case of two scalar fields, though the 
general discussion applies to an arbitrary number of fields. In this case we have 

^ _ C 2 01 + Ci0 2 ^ _ cifa - C 2 4>2 ^ 

'4 + C\ ' y/cf + C 2 

and we can write the explicit dependence of the potential (jSJ) on the orthogonal field [3H [22] 

U(x) = Vi e^ Cl/c2)cx + V 2 e {c2lci)cx . (8) 

This is clearly a positive function (given V± > and V 2 > 0) and bounded from below. Close 
to its minimum we can expand this as [MJ [22] 



U( X ) = U 



c 2 

l + -(x-Xo) 2 + -. 



(9) 



Thus we can confirm that there is a classical trajectory for the two fields in which x 
remains fixed, x — Xo ? while the adiabatic field a rolls down a steep exponential potential 
([S]), which reduces to 

V\ x=X0 = -U e- c % (10) 

and we identify the scaling solution (J3]) with the usual power law solution for a single field 
solution for a negative exponential potential which exists for c 2 > 6 [TU [T5] . 
However we also see that the entropy field, Xi has a negative mass-squared 
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Note this effective mass-squared is a function only of the same parameter c that appears in 
the exponential potential for o and not any other combination of the individual values C\ 
and C2- 

Relative to the Hubble scale we have 

m l 2 9 . . 

* 5 3# = -r +2e - (12) 

where we define the dimensionless parameter 

e = ~Jp> (13) 
and we have e = 1/p for a power-law solution, and e ^> 1 during ekpyrotic collapse. 

Thus there is a strong tachyonic instability associated with the ekpyrotic scaling solution 

0- 

There is a simple explanation for why the scaling solution is unstable. This represents 
evolution along an effective exponential potential §5$) which is less steep than any of the 
individual potentials d2J) from which it was constructed: 

c 2 <c, 2 V!. (14) 

This is what makes assisted inflation solutions stable in an expanding universe as the field 
direction with the flatest potential tends to win out at late times. In a collapsing universe, 
by contrast, the field direction with the steepest (negative) gradient and hence the stiffest 
equation of state tends comes to dominate at late times. 

A full phase space analysis is presented in an Appendix. In addition to the scaling 
solution we have fixed points corresponding to any one of the original fields fa dominating 
the energy density where the other fields have negligible energy density. These correspond 
to power-law solutions where 

p=-|, (15) 

for cf > 6. We find that any of these single dominant field solutions is a stable local attractor 
at late times during collapse. We present the phase-space analysis in Appendix A. 

Finally we note that the scaling solution (j3J) is also not the early time attractor. Restrict- 
ing our attention to solutions along the adiabatic field direction, where the entropy fields 
remain in the minimum value of U(xi), we can use the phase-plane analysis for a single scalar 
field with a steep negative potential [TJ]. These tend to be dominated at early times in a 
collapsing universe by the kinetic energy of the scalar fields and the potential is negligible. 
This corresponds to a power law solution this time with 

P=\- (16) 

This is the original pre big bang proposal [1] for an early collapse phase dominated by the 
kinetic energy of effectively massless moduli fields (see also [TTl 135]). 

In summary, we confirm that there exists a scaling solution corresponding to the combined 
evolution of two or more scalar fields rolling down their exponential potentials. However this 
solution is a saddle point in the classical phase-space of the system (see Fig. 1 in Appendix 
A), representing either a set of measure zero in terms of initial conditions, or a transient 
evolution between the early- and late-time attractors that we have identified. 
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3 Perturbations 



For a collapse model where a < 0, small-scale fluctuations that start far inside the Hubble 
length, H~ l , contract less slowly than the Hubble length and may enter an overdamped 
long- wavelength regime on super-Hubble scales, just as sub- Hubble modes are streched up 
to super-Hubble scales during an inflationary expansion. Thus initial vacuum fluctuations on 
small scales can provide the initial state for large (super-Hubble) scale perturbations during 
collapse. 

We will consider the evolution of isocurvature modes in each of the power-law background 
solutions in the previous section. 

Linear perturbations in fields orthogonal to the background trajectory are decoupled from 
first-order metric perturbations so long as the trajectory remains straight in field-space. They 
obey the wave equation for a massive field in an unperturbed FRW metric 



5 X + 3H5 X + 



'k 2 



+ Wly I 8\ — 



(17) 



Note that such isocurvature perturbations are automatically gauge- invariant. Introducing 
the rescaled field v = aSx, and writing the wave equation in terms of conformal time r = 
/ dt/a, we have 



v" + 



; 2 " i 2 2 

k h m v a 

a x 



v = 0. 



where a prime denotes derivatives with respect to r. For any power- law solution a oc \t\ p we 
have 

aH= , , (19) 



It 



and 



_ = _( e _ 2 )a 2 # 2 . 



(20) 



Thus for any isocurvature field we can write 



k 2 + £+ 3% - 2 



l)2 T 2 



v = 



(21) 



where we have used the dimensionless fast (or slow) roll parameters e = 1/p and rj x , defined 
in Eq. (fT2|). 

Using the usual Bunch-Davies vacuum state to normalise the amplitude of fluctuations 
at early times, we obtain the 



7r e 



l(H-l/2)2 



2 fcV2 

where the order of the Hankel function is given by 

2 _ 9 2e 2 - 3e + 3r? 
V ~4 (e- l) 2 



(22) 



(23) 



5 



and we take v > without loss of generality. At late times, — kr ^> v, this yields [13] 

V Sx ^^\ = C^(-kr)^, (24) 

where C u = 2 u ^ 3 ^ 2 T(u)/7e 3 ^ 2 , and hence we obtain a spectrum of isocurvature perturbations 
on large, super-Hubble, scales with spectral tilt 

An ^^hIF = 3 - 2 - (25 > 

3.1 Scaling solution 

We have seen that the isocurvature field has a tachyonic instability with a negative effective 
mass-squared, given in Eq. (fl2|) . Substituting this into the wave equation (I2T1) we obtain 

Thus from Eq. (I25p we obtain a scale-invariant spectrum for the isocurvature perturbations 
either as e — > 0, corresponding to slow-roll inflation, or as e — ■> oo, corresponding to ekpyrotic 
collapse. To leading order in a fast-roll expansion (e ^> 1) we obtain 

An Sx ~ 2 - . (27) 



Thus we obtain a slightly blue spectrum for a steep exponential potential 
but becoming scale- invariant as e —>■ oo. Any deviations from an exponential potential for 
the adiabatic field, and hence the exact power-law collapse, introduces corrections into the 
tilt [25l [26] . but we note that such corrections could also alter the effective mass rj x and the 
tilt becomes model-dependent. 

3.2 Single-field dominated solution 

As remarked earlier, the scaling solution during collapse ([3]) is unstable with respect to col- 
lapse dominated by any single field, (pi. In this case fluctuations in any orthogonal field, 4>j, 
correspond to isocurvature perturbations decoupled from 0j and first-order metric perturba- 
tions, so we can again use Eq. ffTT|) where we identify 5\ = b4>j- However in this case the 
effective mass of the isocurvature field, 

is much less than the Hubble rate, H 2 . Substituting the dimensionless mass parameter, 
r/ x ~ 0, into Eq. (I2TI) . we find 

v 2 ^\, (29) 

and thus the spectral tilt 

An 2 = 2 . (30) 

We recover the same steep blue spectrum for isocurvature perturbations during (stable) 
ekpyrotic collapse dominated by a single-field, as is found for perturbations in the dominant, 
adiabatic field [7]. 
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3.3 Kinetic-dominated solution 



We have seen that an early time attractor is the kinetic dominated evolution with p = 
1/3 proposed in the pre-big bang scenario. In this limit all the scalar field potentials are 
negligible and all field perturbations, adiabatic and isocurvature, acquire steep blue spectra 
with An = 3, if they have canonical kinetic terms [T7] . 

On the other hand there are pseudo-scalar axion fields with non-minimal kinetic terms. 
Such isocurvature fields with a kinetic coupling to dilaton-type fields can acquire a scale- 
invariant spectrum of perturbations during kinetic dominated collapse as first noted in 
Ref. US] (see also pi]). 

4 Discussion 

We have seen that in this model we are only able to generate a scale-invariant spectrum of 
isocurvature field perturbations during collapse in the scaling solution (j3J) which is a saddle 
point in the phase space of the system [23]. Thus it is neither the early, nor late time 
attractor of the system. 

As a result it does not solve the classical homogeneity problem. Large scale inhomo- 
geneneities in the initial spatial distribution of the fields will grow during the collapse phase, 
in contrast to the original single-field ekpyrotic collapse. 

Even if some preceding phase drives the classical background solution to this unstable 
fixed point throughout space, one might worry that quantum fluctuations could destabilise 
this classical evolution. We see from Eq. that for long wavelength perturbations we 
have 5x oc r l ~ 2u /a 2 and in the ekpyrotic scaling solution where v ~ 3/2 this gives 5\ ~ H. 

Indeed it is easy to see that we always require an instability of this form in a canonical 
scalar field in order to generate a scale invariant spectrum during collapse, as the amplitude 
of field perturbations at Hubble exit are also of order H and thus we require the super-Hubble 
perturbations to grow at the same rate to maintain a scale-invariant spectrum. Precisely 
the same form of instability appears in the power-law collapse model with p = 2/3 which 
produces a scale-invariant spectrum of comoving curvature perturbations from quantum 
fluctuations [TTl H21 ESI H3] ■ 

As H is rapidly growing in the ekpyrotic scenario (e ^> 1) this corresponds to a rapid 
tachyonic growth of the entropy field. On the other hand one should remember that the 
minimum duration of such a phase is also set by the rate of change of H. The change in the 
Hubble rate H determiness the range of comoving scales that exit the Hubble scale (as the 
scale factor a changes only slowly). Thus it might be sufficient for \H\ and thus 5\ to grow 
by a factor of fc max /&min ~ 10 5 while observable scales exit the Hubble scale. 

The actual density of long wavelength quantum fluctuations (per logarithmic range in 
wavelength) either from the kinetic energy or potential energy density of the entropy field 
fluctuations, can be estimated relative to the total energy density to be 

Vt x ~ eH 2 . (31) 

Thus the fractional density of the entropy modes remains small and does not threaten to 
destabilise the classical evolution so long as the Hubble rate remains well below the Planck 
scale (here set to unity). 
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In any case we do not necessarily expect the scaling solution, or any classical solution 
to hold arbitrarily close to the Planck scale, especially if we want to appeal to high-energy 
corrections to produce a bounce in the scale-factor, as envisaged in pre big bang models 
[H [35] and now proposed in ekpyrotic models [26j [27]. 

We note that the existence of a natural turning point in the field-space trajectory due to 
the instability of the scaling solution offers the possibility of converting the scale-invariant 
spectrum of isocurvature field perturbations into a scale-invariant spectrum of curvature 
perturbations before the singularity or bounce without introducing any additional mechansim 
[251 [261 [27]. 

Thus the ekpyrotic collapse with multiple fields offers an alternative scenario for the 
origin of structure in our Universe to contrast with inflationary models. However to do 
so we seem to require at least four distinct phases in the evolution: (1) an initial phase 
to establish spatial homogeneity, (2) a scaling collapse phase to produce a scale-invariant 
spectrum of field perturbations, (3) a phase to turn isocurvature field perturbations into 
curvature perturbations, and (4) a bounce phase to turn collapse into expansion. Note that 
phases (3) and (4) could happen in the reverse order, as in the curvaton mechanism. By 
comparison, the simplest inflationary models require only a single quasi-de Sitter phase to 
establish spatial homogeneity plus a scale-invariant spectrum of curvature perturbations. 



A Appendix 

In this appendix, we present the phase space analysis for mutli-fields. Note that a similar 
analysis was done in Ref. [23J for the case c\ = C2 in an expanding universe. The field 
equations and an acceleration equation are 

4/1 + 2h<P' i + a 2 c^e" c ^ = 0, (32) 

^4en?-«v' a ) 1 (33) 

6 i 

where prime denotes a derivative with respect to a conformal time and h = a' /a. The 
Friedman equation is given by 

h 2 = fy? ~ ^e-^*) . (34) 
We define the dimensionless phase-space variables 

Xt = vk__ (35) 

ay/Vie-**' , . 

m = (36) 

As we are interested in a contracting universe, < 0. The field equations and the evolution 
equation for a give the first order evolution equations for the phase space variables 

^ = -3xi(l -CiJ-Vi, (37) 

3 
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^ = Vi \ZY, X ) - C ^\ X ^J > ( 38 ) 

(39) 

where iV = log a. The Friedmann equation gives a constraint 

(40) 

There are n + 2 fixed points of the system where dxi/dN = dyt/dN = 0. 

A: ES 2 = 1 ' Vi = 0- ( 41 ) 



5i : Xi = -^=, yi = -^j-l, Xj = y j = 0, (for j ^ i), (42) 

R ^ 



3p Cj ' J ^ Cjp \ 3p 



i-lV (43) 



where 



P = E4. (44) 

A is a series of kinetic dominated solutions where the potential energy can be neglected. 
Bi is a single field scaling solution where one of the fields is dominating the energy density. 
B is the scaling solution which is realized by the combined evolution of multiple fields. Note 
that at the fixed point B, the solution for the scale factor is given by 

a = (-t) p = (-r)- p/ ^ p - 1 \ (45) 

To study the stability of the solutions, the linear perturbations around each fixed point 
have to be analyzed. For simplicity we restrict ourselves to the case of two fields. We can 
eliminate y 2 using the constraint y\ = — 1 + x\ + x\ — y\. Then the phase space becomes 
three-dimensional space spanned by (xi, x 2 , yi)- There are four fixed points A, Bi,B 2 and 
B. The solution for the linearized perturbations are characterized by three eigenfunctions: 

5x = ■u 1 exp(m 1 A r ) + u 2 exp(m 2 N) + w 3 exp(m 3 iV). (46) 

Stability during collapse (N — > — oo) requires that the real parts of all the eigenvalues m.j 
are positive. The eigenvalues are given by 

. fa 

(47) 
(48) 
(49) 
(50) 



A : 


mi 


= 0, m 2 = 3 


- ^cixi, m 3 = 6 - V6c 2 x 2 , 


B 1 : 


mi 


= cl, m 2 = - 


2 2 

-3 + |, m 3 = -3 + |, 


B 2 : 


mi 


= y, m 2 = 


r 2 2 

-3 + — , m 3 = -3 + —, 


B : 


mi 


= ->-^ 


, m 2 = m + , m 3 = m_, 
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where 

The kinetic-dominated solutions A are unstable for C\X\ > \/6 and C2X2 > y/E. The single 
field scaling solutions, B^, exist for q > Vd and it is always stable. The scaling solution B 
exists for p < 1/3 and there is always one negative eigenvalue, which means that the scaling 
solution is a saddle point. The instability mode is given by 

5x oc e m ~ N oc (-t) pm - . (52) 

For small p < 1, m_ = —1/p, thus 5x oc (— t)^ 1 oc H. This is exactly the instability that 
the entropy field exhibits around the scaling solution. In fact, in terms of e = 1/p, m_ is 
written as 

m_ = -1(3 - e) + ^/ 3 (3e-l)(e-3). (53) 

Then we can show that the solutions for the entropy field Eq. (j2"2l) is given by 

5 X oce Nm ~. (54) 

Figure 1 shows the phase space trajectories around B. We see that all solutions go to 
the stable fixed points B\ or B2 at late times. It is clear that B is a saddle point and any 
solutions near B go to B\ or B2 at late times. This is the origin of the instability of the 
scaling solution. 
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